ApJ, accepted 

Preprint typeset using I^T^X style cmulatcapj v. 10/09/06 



PROPER MOTIONS OF PSRS B1757-24 AND B1951+32: 
IMPLICATIONS FOR AGES AND ASSOCIATIONS 

B. R. Zeiger 1 , W. F. Brisken 2 , S. Chatterjee 3 ' 4 , W. M. Goss 2 

ApJ, accepted 

ABSTRACT 

Over the last decade, considerable effort has been made to measure the proper motions of the 
pulsars B1757— 24 and B1951+32 in order to establish or refute associations with nearby supernova 
remnants and to understand better the complicated geometries of their surrounding nebulae. We 
present proper motion measurements of both pulsars with the Very Large Array, increasing the time 
baselines of the measurements from 3.9 yr to 6.5 yr and from 12.0 yr to 14.5 yr, respectively, compared 
to previous observations. We confirm the nondetection of proper motion of PSR B1757— 24, and our 
measurement of (// Q , /is) = (— 11±9, — 1±15) mas yr -1 confirms that the association of PSR B1757— 24 
with SNR G5.4— 1.2 is unlikely for the pulsar characteristic age of 15.5 kyr, although an association 
cannot be excluded for a significantly larger age. For PSR B1951+32, we measure a proper motion 
of (iJ, a , — (—28.8 ± 0.9, —14.7 ± 0.9) mas yr -1 , reducing the uncertainty in the proper motion by 
a factor of 2 compared to previous results. After correcting to the local standard of rest, the proper 
motion indicates a kinetic age of ~51 kyr for the pulsar, assuming it was born near the geometric center 
of the supernova remnant. The radio-bright arc of emission along the pulsar proper motion vector 
shows time-variable structure, but moves with the pulsar at an approximately constant separation 
~2.5", lending weight to its interpretation as a shock structure driven by the pulsar. 
Subject headings: pulsars: individual (PSR B1757-24, PSR B1951+32) — ISM: individual (G5.4-1.2, 
CTB 80) — stars: neutron — supernova remnants 



1. INTRODUCTION 

Neutron stars (NSs) are born from the core collapse 
of massive stars, whose deaths are marked by super- 
nova remnants (SNRs). Measuring the proper motions 
of young NSs allows their trajectories to be traced back, 
and can thus provide strong tests of proposed NS-SNR 
associations. In cases where such associations can be 
confirmed, an independent age estimate can be derived 
for both the NS and its associated SNR. 

However, such an exercise is subject to several pitfalls. 
On one hand, a typical radio pulsar might be detectable 
for > 10 7 yr, while the associated SNR fades into the 
interstellar medium (ISM) in < 10 5 yr, leaving behind 
young pulsars with no detectable SNRs. On the other 
hand, there are also several SNRs with no associated pul- 
sars. Not all NSs are radio pulsars, but when young, they 
are all likely to be bright in thermal X-ray emission, as 
illustrated, for exa mple, by the detection of p oint sources 
in the SNRs K esJ3_ jVasisht fc GottheUl997h and Kes 79 
(jSeward et al.l l200"3h . The absence o f such X-ray detec- 
tions in a set of nearby, young SNRs ( Ka plan et al.|[2004l 
2006) aggravates the problem. While unusual cooling 
scenarios may be required for NSs in these young SNRs, 
a partial explanation for the "missing" NSs might lie in 
their birth velocities. Radio pulsars have characteristic 
birth velocities ~ 400-5 00 km s _1 (jArzoumanian et al.l 
l2002t iHobbs et al.ll2005f ). and they are likely to escape 
their natal remnants once the expansion of the SNR is 
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decelerated by the ISM, making identification of pulsar- 
SNR pairs problematic. 

Thus, cases where a young radio pulsar can be as- 
sociated with a SNR are particularly important. The 
pulsars B1757— 24 and B1951+32 present two such sit- 
uations in which associations might be possible. Both 
are young pulsars with high spin-down energy-loss rates, 
and in both cases, the ram pressure balance between 
the NS relativistic wind and the ISM produces a pul- 
sar wind nebula (PWN) with a bow shock structure. 
PSR B1757— 24 appears to be exiting the approximately 
circular SNR G5.4— 1.2, whose asymmetric brightness, in 
conjunction with the PWN produced by the pulsar, pro- 
duces the structure known as "the Duck" (see Fig. QJ. 
Meanwhile, PSR B1951+32 drives a complex interaction 
within the SNR CTB 80 (G69.0+2.7; see Fig.©, includ- 
ing a radio-bright structure resembling a bow shock in 
the direction of its motion. 

We present Very Large Array (VLA) observations that 
improve the time baseline for the proper motion mea- 
surements of these two pulsars. In Sj2]we detail the ob- 
servations, analysis, and results for PSR B1757— 24, and 
likewise for PSR B1951+32 in £3 We discuss the impli- 
cations of our results in ijll 

2. PSR B1757-24: OBSERVATIONS AND RESULTS 

PSR B1757— 24 is an energetic young pulsar, with 
a period P = 125 ms, a spindown energy loss rate 
E = 10 36 ' 4 erg s _1 , and a characteristic age r c = 
P/2P = 15.5 kyr. As shown in Figure [H its location 
and the morphology of its PWN suggest that it is escap- 
ing the circular SNR G 5.4-1.2 (|Frail fc Kulkarnil Il99lt 
iManchester et al] fl99lK If it was born at the center of 
the remnant and its characteristic age is close to its real 
age, that would imply a proper motion ~ 70 mas yr -1 
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and a transverse velocity —1500-2000 km s , high 
compared to the radio pulsar population. However, 
iGaensler fc Fraill (|2000t ) placed an upper limit on the 
motion of the PWN that was inconsistent with such 
a high velocity, and suggested that the true age of 
the pulsar was much larger than the characteristic age. 
iThorsett et al.l ()2002f ) derived comparable limits on the 
proper motion 5 of the pulsar (rather than the PWN) of 
(fi a ,fis) = (— 2.1±7.0, — 14±13) mas yr -1 and suggested 
instead that the proximity of B1757— 24 and G5.4— 1.2 
was merely a line-of-sight coincidence in the crowded re- 
gion near the Galactic center. Arguments for and against 
an associati on between the pul sar and SNR were sum- 
marized by iBlazek et al.l ((2006), who improved the up- 
per limit on the westward proper motion of the PWN 
and concluded that both interpretations for the Duck (a 
line-of-sight coincidence and a genuine association with 
a large age) remained viable. 

Based on the dispersion measure, PSR B1757— 24 
is estimated to lie at a distance of 5.2±0.5 kpc 
(jCordes fc Lazioll2002D . while the distance to G5.27-1.2 
is greater than 4.3 kpc based on H I absorption 
(jFrail et al.lll994t IThorsett et~alll2002f l. Here we adopt 
a distance of 5d§ kpc. 

2.1. Observations 

VLA observations of PSR B1757-24 spanned 6.4 yr, 
from 1998 to 2004. The 1998 June 03 epoch was observed 
using the BnA configuration; all of the later epochs used 
the A configuration. At the frequency of observation, 
1.45 GHz, the A configuration resolution at a declina- 
tion of -25° is -1.2" x 2.2", elongated north-south. The 
first three epochs employed the pulsar gate on the right 
circular polarization to allow increased sensitivity in the 
pulsar measurements by only accepting data when the 
pulsar was "on." The 2002 February 21, 2002 May 02, 
and 2004 December 09 epochs included the Pie Town 
VLBA antenna, more than doubling the maximum east- 
west baseline length. The ungated data from the 2002 
February 21 epoch was not used, as the observation was 
not long enough to provide sufficient u-v coverage and 
sensitivity needed for clear identification of the pulsar. 
In order to minimize vari ation in u-v coverage between 
epochs and in contrast to IThorsett et al.l (|2002f) . we do 
not use the Pie Town baselines. Observations were made 
in the '2AD' correlator mode, offering 15 spectral chan- 
nels each of width 1.56 MHz in both circular polariza- 
tions and allowing a usable field of view of —30'. Details 
of these observations are shown in Table Q] 

The data were calibrated using standard procedures 
within the AIPS 6 package. Sources 3C 286 and 3C 48 
were used as flux density and bandpass calibrators and 
1751—253 was used for phase calibration. The AIPS task 
UVFIX was used to recalculate the UVW coordinates in- 
corporating aberration corrections not performed by the 
VLA online system. Given the large field of view and the 
sparse distribution of sources with compact structure, we 
selected 14 sub-regions of the field (containing the pulsar 
and 13 other apparently compact sources) for imaging. 
These fields were jointly deconvolved to produce initial 

5 All uncertainties are 68% confidence intervals, except where 
both 68% and 95% intervals are explicitly reported. 

6 See http://www.aips.nrao.edu/ 



images. Two self-calibration and imaging iterations were 
performed using all 14 sub-regions to generate the final 
images used for the astrometry. 

Positions of reference sources and the pulsar were de- 
termined with a Gaussian fit using the AIPS task JM- 
FIT. The proper motion of the pulsar was measured with 
respect to a reference frame defined by six point sources 
chosen based on compactness from the 13 imaged sub- 
regions. The reference source positions and fluxes from 
the (ungated) 2004 December 09 epoch are listed in Ta- 
ble [2j and our met hod is described in greater detail in 
iMcOarv et all (l2?ffllh . 

The 10" PWN in which the pulsar is embedded is re- 
solved out by our observations, but it contributes an ad- 
ditional uncertainty to the position of the pulsar in the 
ungated data. To account for the ~50% increase in rms 
noise within this region relative to the background, the 
pulsar position uncertainties were increased in ungated 
epochs by a factor of 1.5. In the gated data, the PWN 
contribution to the noise is insignificant and the posi- 
tion uncertainties of the reference sources dominate the 
proper motion uncertainty. 

2.2. Results 

The structure of the PWN around the pulsar is un- 
certain and could cause a systematic offset in the mea- 
surement of the pulsar position, affecting registration be- 
tween the gated and ungated data sets. Thus, we mea- 
sure the proper motion of the pulsar separately for the 
gated and ungated data and combine the two indepen- 
dent results for a final proper motion value. Using un- 
gated data, we measure a proper motion of (fx a ,fj,s) = 
(-16.9 ± 12.6, -35.8 ± 22.4) mas yr" 1 ; with gated data, 
we measure ([i a ,[is) = (— 5.9±11.8, 6.0±20.3) mas yr -1 . 
Combining the results gives a measurement of (/x Q , ^s) — 
(—11 ±9, — 1 ± 15) mas yr -1 , consistent with a nondetec- 
tion. Note that correction to the local standard of rest 
(LSR) for the pulsar involves adjustments at the level of 
1 mas yr . Since these are much lower than the uncer- 
tainties, we ignore these corrections for B1757— 24. 

We note that our results and those of IThorsett et all 
(2002) arc not independent, since two epochs of data 
are common between the two analyses, although the Pie 
Town link data were not used i n the present work. The 
measurement of IBlazek et al.1 (|2006l ) is completely in- 
dependent, since they observe the nebula and not the 
pulsar. However, only one dimension of motion was 
probed by these authors. We derive a 68% confidence 
upper limit on the westward motion of the pulsar, [i a > 
— 14.9 mas yr -1 (fi a > —24.8 mas yr -1 at 95% confi- 
dence), which corresponds to a westward transverse ve- 
locity v± < 360d 5 km s -1 (v ± < 600d 5 km s _1 at 95% 
confidence). The implications are discussed further be- 
low, but we note that such velocity limits are consistent 
with th e velocity distribution of ordinary young ra dio 
pulsars (|Arzoumanian et al.ll2002l ; iHobbs et alJl2005T ). 

3. PSR B1951+32: OBSERVATIONS AND RESULTS 
PSR B1951+32, a 39.5 ms pulsar with r c = 107 kyr 
and a spin-down energy-loss rate E — 10 36 - 6 erg s -1 , lies 
near the southwestern edge of the approximately circu- 
lar infrared shell of CTB 80 (jFesen et al.l fl988l ) . The 
pulsar is surrounded by a ~30" diameter asymmetric 
PWN at the western edge of an 8'x4' east- west plateau 
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Fig. 1. — The Duck: the SNR G5.4-1.2, G5. 27-0. 90, and B17 57-24 system. Lef t: Large-scale structure of the remnant, as observed 
by the VLA at 327 MHz with an ang ular resolution of 40" X 23" (Brogan ct al. 2006, image data courtesy of the authors). The extended 
circular structure is SNR G5.4— 1.2, while G5. 27— 0.90 is the "head" of the Duck protruding from the west edge of the remnant. The 
small white square shows the approximate location of the right-hand panel. Right. Pulsar B1757— 24 and its PWN (the western tip of 
G5. 27— 0.90), as observed at 1.4 GHz in 2002. The large-scale structures are resolved out by the interferometer response, leaving only the 
"beak" visible at an angular resolution of 2" 2 X l'/2. 



TABLE 1 

Observational Parameters for B1757— 24 



Epoch 


Obs. 


Poln. 


Gated 


Flux density 


rms noise 


Beam 


Pos. 


Tint 




Code 






(mjy) 


(mjy beam ) 


(") 


angle 


(hr) 


1998 Jun 03 


AF336 


R 


Y 


5.0 


0.3 


2.3x1.2 


16° 


1.7 






L 


N 


1.9 


0.1 


2.2x1.2 


16° 


1.7 


2001 Jan 01 


AB969 


R 


Y 


3.1 


0.2 


2.3x1.2 


-1° 


3.2 






L 


N 


1.1 


0.1 


2.2x1.2 


0° 


3.2 


2002 Feb 21 


AB1029 


P, 


Y 


3.4 


0.2 


2.0x1.2 


5° 


1.5 


2002 May 02 


AB1029 


1 


N 


4.5 


0.1 


2.2x1.2 


4° 


2.3 


2004 Dec 09 


AB1139 


I 


N 


1.3 


0.1 


2.3x1.2 


0° 


3.5 



Note. — All observations are at a frequency of 1.45 GHz. Pulsar gating improves the apparent 
pulsar flux density by a factor < [(T on + T^/Ton] 1 / 2 . 



surements, while IMigliazzo e t al. (2002) directly measure 
a proper motion /i = 25 ± 4 mas yr _1 at a position angle 
252° ± 7° after correcting for the effects of differential 
Galactic rotation. Their implied transverse velocity of 
the pulsar is V± = (240 ± 40) d 2 km s" 1 . 

I n the direction o f the p ro per motion vector m easured 
by IMigliazzo et~aH (|2002D . iMoon et alj (|2004D observe 
a cometary X-ray synchrotron nebula and an Ha bow 
shock. The X-ray emission peaks at the pulsa r loca- 
tion and is confined within the Ha structure (|Hesterl 
2000) that is clearly defined at an angular separation of 
~7" from the pulsar. The X-rays are produced by syn- 
chrotron emission from the pulsar wind, confined within 
a bow shock produced by the ram pressure of the wind 
interacting with the SNR wall. 

In radio images, B1951+32 appears to be situated just 
inside a limb-brightened bubble ~30" in diameter (see 
Fig. [21 right). A 5" portion of the bubble nearest to the 
pulsar is substantially brighter than any other portion 
of the shell. In MERLIN observatio ns at 1.6 GHz with 
0.15" resolution, Golde rTet alj ([2005) find that the radio- 
bright arc shows compact structure, resembling a radio 
bow shock. The structure is ~2.5" from the pulsar, which 
puts it at the head of the cometary X-ray nebula, but 
enclosed within the Ha bow shoc k which is ~7" away 
from the pulsar (|Moon et al.ll2004D . 



TABLE 2 

Reference Sources Used in the Proper 
Motion Fit of B1757-24 



R.A. 


Deck 


Flux density 


$scp 


(J2000) 


(J2000) 


(mjy beam ) 


(') 


18 01 27.40 


-25 07 38.1 


8.7(2) 


17.3 


18 01 58.57 


-24 55 48.5 


6.6(2) 


14.0 


18 00 47.31 


-24 42 45.0 


2.6(1) 


9.2 


18 00 15.71 


-25 00 22.0 


3.6(1) 


13.4 


18 01 28.11 


-24 34 23.7 


1.2(2) 


18.2 


18 00 41.14 


-24 42 04.4 


4.9(1) 


10.3 



Note. — Units of right ascension are hours, min- 
utes, and seconds, and units of declination arc de- 
grees, arcminutes, and arcseconds. 



of emission (|Castelletti et al.ll2003l) . CTB 80, observed 
at 1.4 GHz to have three arms that cover 1.8 d eg 2 and 
converge near B1951+32 (|Castelletti et al.|[2003| ). has an 
expanding H I shell that yields a dyna mical age estimat e 
of 77 kyr for the pulsar-SNR system (jKoo et alJli990ft . 
From its dispersion measur e, the distance to B195 1+32 is 
estimated as 3.1 ± 0.2 kpc (|Cordes fc Laz io 2002), while 
the distanc e to CTB 80 is estimated as 2 kpc from H I 
absorption (|Strom fc Stappe rs 2000). Her e we adopt a 
distance of 2d? kpc. iKulkarni et alj (|1988f ) estimated a 
speed of 300 km s" 1 for the pulsar from scintillation mea- 
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Fig. 2.— CTB 80 (G69.0+2.7) and the PWN produced by PSR B1951+32. Left: Infrared ratio image (60^m/100/im) constructed from 
the IRAS archive shows the nearly complete shell of CTB 80. Overplotted lines indicate the projected path of the pulsar, projected back in 
time with the 1 a uncertainty on the measured proper motion, corrected to the LSR. Circles indicate the estimated birth locations for age 
estimates of t p ~ 51 kyr (the age of closest approach to the SNR geometric center, as shown with the white cross) and the ch aracteristic 
(spin-down) age t c = 107 kyr. The projected path does not pass close to the expansion center obtained by IKoo et al.l 1 119901) . indicated 
by the white triangle. Right: VLA 1.4 GHz image of B1951 + 32 and its PWN, with an a ngula r resolution of l'.'0x0'.'9 and rms noise 
0.03 mjy beam -1 . The arrow represents 300 yr of travel along the proper motion vector (see t|3.2j l. The bow shock is visible as the bright 
arc to the south-west of the pulsar. Note the difference in scale between the two panels. 



3.1. Observations 

Measurements of the proper motion of PSR B1951+32 
were made using five epochs of archival VLA data from 
1989.04 to 2003.50. Observation details are shown in Ta- 
ble [31 All observations were made with the VLA A array 
except the 2000 November 25 epoch, which included the 
Pie Town antenna. To maintain consistent u-v coverage 
in each epoch, the baselines to Pie Town were removed 
before data were reduced. Observations were made with 
15 channels of 1.56 MHz each. Data were calibrated 
with VLA calibrator sources 3C 286 (for flux density 
and bandpass calibration) and J1925+211 (for phase cal- 
ibration), following the method outlined in H2.1[ Nine- 
teen sources with compact structure were found in the 
field of view and were included in the imaging and self- 
calibration process. Eight pointlike sources (see Table |3J 
were used in the determination of the proper motion in 
the manner described in §2.11 

The pulsar B1951+32 lies in a complicated region of 
the CTB 80 SNR and is embedded in its PWN, making 
the position fit sensitive to observation parameters such 
as observing frequency and u-v coverage at each epoch. 
To account for this to first order, a linear brightness gra- 
dient in a 1.75" square around the pulsar was fit simul- 
taneously with the pulsar position. Trial fits performed 
without this procedure resulted in substantially greater 
scatter in the position measurements. Additional trials 
demonstrated that the fit was not very sensitive to the 
exact size and position of the region fitted. 

3.2. Results 

The proper motion of B 195 1+32 was measured as 
= (-28.8 ±0.9, -14.7 ±0.9) mas yr" 1 . Figured 
demonstrates the significance of the proper motion rel- 
ative to that of the eight presumably stationary sources 
used to define the frame in which the proper motion was 
measured. For each of the reference sources, we derive 
a proper motion in the frame defined by the seven other 
sources. As shown in Figure [3j these proper motion val- 
ues are consistent with zero at ~1 a, indicating that a 
self-consistent reference frame has been established. 
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Fig. 3. — Proper motion of PSR B1951+32 and the reference 
sources used for the fit. Reference source motions with respect to 
each other are consistent with noise with and scattered around zero. 
The pulsar lies in the bottom right corner with a proper motion of 
(H„,/1 S ) = (-28.8 ± 0.9, -14.7 ± 0.9) mas yr -1 . In this figure the 
derived motion of the pulsar is relative to the frame defined by all 
eight reference sources. 

A significant correction is required to obtain the mo- 
tion of the pulsar in its LSR so that it can be related to 
the remnant. We use a flat Galactic rotation curve with 
velocity 220 km s _1 and assume a solar radius of 8.5 kpc. 
Further correction for the solar peculiar motion is applied 
(|Dehnen fc Binnevlll998 l ) . The magnitude of the correc- 
tion depends on the distance; an additional 0.6 mas yr -1 
has been added in quadrature to the proper motion 
uncertainty to accommodate the magnitude range of 
the correction over the 1-3 kpc distance range consid- 
ered. The LSR-corrected proper motion is found to be 
(a4X) = (-26.9±l.l,-10.5±l.l)mas yr _1 , for a total 
proper motion \i! = 29 ± 1 mas yr -1 at a position angle 
249 ± 2° east of north. 

3.3. Shock Separation 

To maximize sensitivity to the extended structure of 
the shock, all u-v data within each epoch were added with 
the AIPS task DBCON. Before reimaging, the resulting 
data were convolved with a circular beam with a 1.2" 
diameter, the longest semi-major axis of any beam in any 
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TABLE 3 

Observational Parameters for B1951+32 



Epoch 


Obs. 


IF 


Pol. 


Froq. 


Flux density 


rms noise 


Beam 


Position 


T int 




Code 






(MHz) 


/ T \ 

(mjy) 


(mjy bm 1 ) 


( ) 


Angle 


(hr) 


1989 Jan 13 


AS357 


1 


RR 


1385 


1.01 


0.07 


l.lxl.l 


50° 


5.9 






2 


LL 


1652 


0.62 


0.09 


1.0x0.9 


81° 


5.9 


1991 Jul 18 


AF214 


1 


RR 


1385 


1.95 


0.08 


1.3x1.2 


-47° 


5.8 






2 


LL 


1652 


1.32 


0.08 


1.0x0.9 


-74° 


5.8 


1993 Jan 08 


AF235 


1 


RR 


1385 


1.75 


0.06 


1.2x1.2 


-64° 


5.9 






2 


LL 


1652 


1.79 


0.05 


1.1x1.0 


-58° 


5.9 


2000 Nov 25 


AG602 


1 


RR,LL 


1385 


0.63 


0.05 


1.3x1.2 


-69° 


5.2 






2 


RR,LL 


1516 


0.91 


0.07 


1.1x1.0 


-75° 


5.2 


2003 Jun 03 


AG650 


1 


RR,LL 


1665 


1.05 


0.03 


1.0x0.9 


-70° 


6.2 



TABLE 4 

Reference Sources Used in the Proper 
Motion Fit of B1951+32 



R.A. 


Decl. 


Flux density 


$scp 


(J2000) 


(J2000) 


(mjy beam ) 


(') 


19 53 23.70 


+32 53 35.3 


4.7(3) 


5.4 


19 53 16.37 


+32 48 46.5 


21.7(3) 


5.5 


19 53 17.92 


+32 49 11.4 


6.8(3) 


5.4 


19 53 13.42 


+33 01 21.8 


97.8(3) 


9.3 


19 52 45.25 


+33 03 30.0 


6.2(3) 


11.2 


19 52 15.80 


+32 49 35.7 


839.3(8) 


9.4 


19 52 30.01 


+32 55 26.2 


5.4(3) 


6.5 


19 52 13.32 


+32 59 22.1 


61.7(3) 


11.6 



Note. — Units of right ascension are hours, min- 
utes, and seconds, and units of declination are de- 
grees, arcminutes, and arcseconds. 



of the concatenated data sets, to ensure equal sensitivity 
to large-scale structure across all epochs. 

Full hydrodynamical modeling of the shock is beyond 
the scope of the present work. Instead, we simply ex- 
tract the flux density in slices along the observed proper 
motion vector in the LSR (position angle = 249° east of 
north) in order to estimate the angular separation be- 
tween the pulsar and the shock structure at each epoch. 
As shown in Figure [H the pulsar position is well defined 
in these slices, but the transverse structure of the shock 
changes with epoch and the separation of the pulsar from 
the shock peak (A6*p ca k) varies between 2.4" and 3.1". 
At each epoch, given the map rms noise <r, we define 
the shock thickness as the range where the emission is 
within 1 a of the shock peak. The near and far edges of 
the shock, defined in this manner, are listed in Table [5] 
and indicated in Figured! 

Over the 14.5 yr span of observations, the pulsar moves 
by ~ 0.47", but the distance between the pulsar and 
the near edge of the shock is ~2.1"±0.2" and the sep- 
aration to the far edge is ~3.6"±0.6". Since the sep- 
aration between the pulsar and the radio-bright fea- 
ture does not show a secular decreasing trend, we con- 
cur with previous inferenc es (jHester &; Kulkarnl 119881 : 
IChatterjee fc Corded 120021 ) that the feature must be a 
shock driven by the pulsar wind as the pulsar travels 
through the ambient medium. 

4. DISCUSSION 

4.1. The Association between SNR G5.4-—1.2 and 
PSR B 1757- 24 

The identity of the remnant associated by birth with 
the pulsar B1757— 24 is still a mystery. Based on mor- 
phology and location, SNR G5.4— 1.2 is the only known 




-2 2 
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Fig. 4.— Slices through PSR B1951+32 and its bow shock. The 
flux density is plotted as a function of the offset from the derived 
pulsar position (as indicated by the solid line), along the measured 
proper motion vector. Successive epochs are vertically displaced by 
1 mjy (except for 2003, which is displaced by 1.25 mjy) for clarity, 
but note that the time interval between successive observations is 
not uniform. The map rms noise a at each epoch is indicated by 
the error bar at the far left. The peak of the shock is marked, and 
the width of the shock within 1 a of the shock peak is also shown 
by the horizontal error bar. The measured separation is consistent 
with being constant over ~14.4 yr, while the pulsar moves ~0.47"in 
that time period, as shown by the horizontal arrow. 

remnant likely to be related to the pulsar, but that asso- 
ciation has serious problems (see Blazek et al. [2006] for 
a recent discussion). To remain tenable, such an associa- 
tion requires a pulsar age much larger than the character- 
istic age r c = 15.5 kyr, an explosion site much closer to 
the current location of the pulsar, or some combination of 
the two. For example, G5.4— 1.2 coul d be expanding int o 
a region with a density gradient (e.g.. lGvaramadzef 2004). 
The brightening of the remnant on the western limb, the 
side nearer to the Galactic plane, may be evidence for 
such an asymmetric expansion. However, even in an op- 
timistic scenario, the birth site could be no closer than 
the westernmost limb of G5.4— 1.2. In such a case, for a 
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TABLE 5 
Angular Separation between 
PSR B1951+32 and the Radio 
Bow Shock 



Epoch 


A^Pcak 


A0N 8ar 


A0 Far 




(") 


(") 


(") 


1989 


2.4 


1.9 


3.2 


1991 


3.1 


2.0 


3.5 


1993 


3.1 


2.4 


3.6 


2000 


2.5 


1.9 


4.3 


2003 


2.6 


2.2 


3.0 



Note. — We list the angular dis- 
tance from the pulsar peak to the peak 
of the shock and the range within 1 a 
of the shock peak at each epoch. The 
position fit uncertainties are insignifi- 
cant compared to the thickness of the 
shock itself. 

true age < r c , a westward proper motion > 18 mas yr _1 
is required for an association, but all three measurements 
reported in Table [6l where we summarize our results, im- 
ply westward motions no greater than 15 mas yr" 1 . It 
thus appears that either the true age is larger than r c , 
probably significantly so, or the association of the pulsar 
with G5.4— 1.2 must be excluded. 

We note that the characteristic ages of some 
young pulsars with well-determined kinet ic ages are 
overes timates, for ex ample, J181 1 — 1 925 (jKaspi et al.l 
l2001bH. J0538+ 2817 (iNg et all l200l. and B1951+32 
(Migl iazzo et alJl2002L our $L2]) . Thus, the true age of 
B1757— 24 might be no greater than 15.5 kyr, render- 
ing the association unlikely. But Vela, for example , may 
have a true a ge gre ater than t c (jLvne et al.lll996l ). and 
iBlazek et all |2006) argue that a growing surface mag- 
netic field, which decouples the true age from the char- 
acteristic age, could account for such an age discrepancy 
in the Duck while preserving the association. Given a 
long enough time baseline, the pulsar proper motion (or 
the motion of the nebula) will surely be detectable, al- 
lowing a firm conclusion about this vexing association. 

4.2. The Age and Bow Shock of B 195 1+32 

For B1951+32, we measure a proper motion of 
(Mct.M*) = (-28.8 ± 0.9,-14.7 ± 0.9) mas yr" 1 , con- 
firmin g and improving o n the r esults of iMigliazzo et all 
(|2002ft . IMigliazzo etail (|2002ft filter the u-v plane to 
remove all spatial scales greater than 4" and to leave 
only compact sources. We use a different approach to re- 
move the complexity of the PWN, subtracting the linear 
slope of the PWN from the region around the pulsar in 
the image domain. The inclusion of a longer time base- 
line and a larger field of reference sources increases both 
the pre cision and the ac c uracy of this measurement over 
that of IMigliazzo et all (|2002l ). The new measurement 
corresponds to an LSR-corrected transverse velocity of 
(274 ± 12) d 2 km s -1 . At a position angle of 249°, the 
proper motion vector is well matched to the lon g sym- 
metry axis of the X-ray PWN (jMoon et alJl2004f ). 

The measurement of the proper motion allows a revised 
estimate for the age of the pulsar. The LSR-corrected 
proper motion of the pulsar is shown in Figure [2 In 
projection, it s closest approac h to the expansion center 
estimated by iKoo et ail |l990D is ~330", a 5.5 a devia- 



tion assuming a 1' uncertainty in the reported expansion 
center. A geometric center at coordinates 19h54m50s, 
33°00'30" (J2000.0) was estimated by superposing a cir- 
cle on the remnant shell seen in the left panel of Figure [2j 
The pulsar's closest approach to this point implies an age 
t p ~ 51 kyr with a minimum separation of 74", consistent 
with the 56" positional uncertainty based solely on the 
proper motion uncertainty. Our measurement of the pul- 
sar motion is therefore not consistent with t he ~77 kyr 
age a nd the expansion center estimated by IKoo et aTl 
(1990) for the CTB 80 SN R but it is consis t ent w ith the 
pulsar age determined by IMigliazzo et alJ (|2002f ) . We 
note that the true explosion center may be masked by 
asymmetric expansion due, for example, to density gra- 
dients and structure in the ISM, and neither estimate for 
the age may be correct. However, the true age is almost 
certainly less than r c ~ 107 kyr, since the pulsar reaches 
the far edge of the SNR when its measured proper mo- 
tion is projected backwards over that time interval, as 
shown in Figure [2] 

As B1951+32 moves through its ambient medium, its 
relativistic wind outflow is confined by ram pressure, and 
as seen in some other PWNe, synchrotron emission from 
the confined pulsar wind produces a cometary X-ray neb- 
ula. On the other side of the contact discontinuity, the 
shocked IS M produces Ha emis sion (see, e.g., the illus- 
tration bv lGaensler et al.ll2004l ). In such a picture, the 
radio-bright arc would be produced by the confined pul- 
sar wind at the location of the contact discontinuity. As 
we have shown here, the structure of the bright arc is 
dynamic, changing from epoch to epoch, but overall, it 
moves with the pulsar at an approximately steady sep- 
aration ~2.5", corresponding to a standoff distance of 
~0.024 d 2 pc. 

4.3. Two Similar Pulsars in Very Different 
Environments 

Currently, over 1600 radio pulsars have measured val- 
ues f or P and P, as listed in the ATNF pulsar cata- 
log 7 (Manches ter et aLll2005h . In this large ensemble, 
B1757-24 and B1951+32 are among the fifty most en- 
ergetic pulsars, with E = 10 364 and 10 36 6 erg s _1 , re- 
spectively. They are also among the hundred youngest 
known pulsars, with characteristic ages r c = 15.5 and 
107 kyr, respectively. On one hand, the characteristic 
age of B1951+32 is likely to be an overestimate, since its 
current spin period is not muc h greater than the typical 
birth spin periods of pulsars (jFaucher-Giguere fc Kaspil 
2006), and its kinematic age t p ~ 51 kyr, assuming birth 
near the center of CTB 80. On the other hand, the up- 
per limits on its westward proper motion suggest that 
B1757— 24 may be substantially older than its character- 
istic age. Although the true ages of the pulsars differ in 
relation to their respective characteristic ages, these two 
pulsars are quite similar in energetics and age. 

Both pulsars also exhibit cometary PWNe. B1757— 24 
is associated with a ~ 10" elongated nebula visible 
at both radio (as in Fig. [T]) and X-ray wavelengths 
dKaspi et al.ll2001aft. B 1951+32 shows a cometary PWN 
in X-rays ( Moon et al.ll2004h . but its radio structure re- 
sembles a bow shock and an extended, limb-brightened 
bubble, and it is enclosed by a complex nebula visible 

7 See http://www.atnf.csiro.au/research/pulsar/psrcat 



Proper Motions of PSRs B1757-24 and B1951+32 
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TABLE 6 

Summary of Proper Motion Measurements 



This work 



B1757-24 
Ref. 1 



Ref. 2 



B1951+32 
This work Ref. 3 



Target 


B1757-24 


B1757-24 


G5. 27-0.9 


B1951+32 B1951+32 


y, a (mas yr _1 ) 


-11 ± 9 


-2.1 ± 7.0 


> -7.9 


-28.8 ±0.9 -29 ±2 


AM (mas yr" 1 ) 


-1 ± 15 


-14 ± 13 




-14.7 ±0.9 -8.7 ±1.3 


Time Baseline (yr) 


6.5 


3.9 


12.0 


14.5 11.9 


Frequency (GHz) 


1.1 


1.4 


8.5 


1.4 1.4 


RA of Pulsar 


18 01 00.016±0.008 




19 52 58.206±0.001 


DEC of Pulsar 


-24 51 27.5±0.2 




32 52 40.51±0.01 


Epoch of RA, DEC 


2004 Dec 09 




2003 Jun 03 



References. — (1) IThorsett et~aTI [|2002h ; (2) IBlazek et all ||2U06D ; (3) 

IMigliazzo et~aTI (2002) 

Note. — - All p rovided uncertainties and limits are are 68% confidence intervals. 
IBlazek et all (f200a) do not measure i ig, an d their 68% fi a has been inferred from their 
5 <t limit; results of Migliazzo et al. (2002) have been adjusted to remove the differen- 
tial Galactic rotation applied in the published result. Proper motions and positions are 
reported in J2000.0 coordinates and do not include LSR corrections. Units of right ascen- 
sion are hours, minutes, and seconds, and units of declination are degrees, arcminutes, 
and arcseconds. 



in Ha (jHesterl I2OOO0. with extended lobe s and a bow 
shock structure. iHester fc Kulkarnll (|l988l ) suggest that 



the pulsar is interacting with the wall of its evolved SNR, 
and we find that the radio shock structure changes with 
time while moving with the pulsar. On the other hand, 
Bf757— 24 appears to be outside any parent remnant. 
If indeed it is associated with G5.4— 1.2, it would have 
interacted with the wall of the remnant in the past (as 
B1951±32 is at present), possibly re-energising the shell 
and leaving behind the structure we see as the Duck. Al- 
ternatively, the parent remnant of B1757— 24 is no longer 
visible, and we are witnesses to a cosmic coincidence. 
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